TTP-04-08 
SFB/CPP-04-16 
May 2004 
hep-ph/0405193 



O- Four- loop renormalization of QCD: full set of 

renormalization constants and anomalous 
^ : dimensions 

m : K.G. Chetyrkin 1 

Institut fur Theoretische Teilchenphysik, Universitat Karlsruhe, D-76128 Karlsruhe, 

^ , 

CO . Germany 
G\ 

in 
o . 

■^j- . Abstract 

O , 

The anomalous dimensions of the gluon and ghost fields as well as 
those of the ghost-ghost-gluon and quark-quark-gluon vertexes are 
^1 analytically computed at four loops in pQCD. Taken together with 

already available anomalous dimensions of the coupling constant, the 
quark field and the mass the results lead to complete knowledge of 
^ ' all renormalization constant entering into the renormalization of the 

QCD Lagrangian at the four- loop level. As a by-product we get scale 
and scheme invariant gluon and ghost propagators at NNNLO. Using 
a theorem due to Dudal, Verschelde and Sorella, we also construct the 
four- loop anomalous dimension of the "gluon mass operator", A 2 , in 
the Landau gauge. 
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1 Introduction 



The renormalization group equation is a powerful tool in investigating the proper- 
ties of the Green functions of a renormalizable field theory. Its crucial ingredients 
are the anomalous dimensions of quantum fields as well as those of mass and 
coupling constant (s). 

In recent years there has been achieved a significant progress in perturbative 
calculation of higher orders corrections to renormalization group functions For 
example, the most physically important RG functions of QCD — the (3 function 
and the quark mass and field anomalous dimensions — have been computed at a 
record-setting four-loup level [1-4]. 

In the same time the anomalous dimensions of the gluon and ghost fields are 
available in literature only at three- loop level [4,5]. It is unfortunate for at least 
two reasons. First, the gluon and ghost field anomalous dimensions are important 
in comparison of the non-perturbative results for the momentum dependence of 
the corresponding propagators with perturbative predictions [6-17]. Second, only 
the knowledge of anomalous dimensions of all fields of the QCD Lagrangian leads, 
together with the f3- function, to complete reconstruction of all Renormalization 
Constants (RCs) entering into the renormalization of the QCD Lagrangian (see 
below) . 

In the present paper we fill the gap by analytically computing the anomalous 
dimensions of the gluon and ghost fields as well as that of the ghost-ghost-gluon 
vertex at four loops. All calculations have been done in the general covariant 
gauge. 

We apply our results to find the scheme and scale invariant gluon and ghost 
propagators at Next-Next-Next-Leading Order (NNNLO) as well as the four-loop 
anomalous dimension of the composite operator A 2 . 



2 Notations and generalities 

The QCD Lagrangian with rif quark flavors in the covariant gauge reads: 



£ = ~G^C^ + ^(ip-m f )^ (1) 
4 /=i 

^(PAff + Pridf-gf^JAl), (2) 

where 

G% = d,Al-d u Al + g{A^A u )\ (A x B) a = f abc A b B c , (3) 

D, = d,-igA*T a , p = YD,. (4) 
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The quark field ip- has a mass m/ and transforms as the fundamental repre- 
sentation and the gluon fields as the adjoint representation of the gauge group 
SU(3). and f abc are the generators of the fundamental and adjoint represen- 
tation of the corresponding Lie algebra. The c a are the ghost fields and £^ is the 
gauge parameter (£x, = corresponds to the Landau gauge). 

By adding to (1) all counterterms necessary to remove UV divergences from 
Green functions, one arrives at the bare QCD Lagrangian written in terms of the 
renormalized fields 2 : 

£o = ~Z 3 (d fl A u - d u A,) 2 - l -g Z\ 9 (d,A a u - cU£) (A, x A v f 

-g 2 Z{ 9 {A, x A v f - Zf ^-{d u A,f + Z 3 C d u c (d u c) (5) 

+ g Z? 9 d»c (Axc) + Z 2 J2 + gZf 9 Z 2 l A - Z m m f ) . 

/=i 

Here Z| 9 is expressed through the RC of the gauge fixing parameter £ L as follows 

£l,o = Z^ L , Z 2g = Z 3 /Z^. 

Z 3 , Z 2 , Z% are the wave-function RCs appearing in the relations between the renor- 
malized and bare gluon, quark and ghosts fields, viz. 

A? = y/z~ 3 A«*, 4 = Jz~4, c a = Jz§c a . (6) 

The full set of the vertex RCs 

ZY, V E {3g, 4g, ccg, ^g} (7) 

serve to renormalize 3-gluon, 4-gluon, ghost-ghost-gluon, quark-quark-gluon ver- 
tex functions respectively. 

The Slavnov- Taylor identities allows one to express all four vertex RCs in 

terms of an independent one, Z„ = — , and the above listed wave function 
RCs. The corresponding relations are: 



Z ( = Z 3 , (8) 

z 9 = s[zj 9 {z 3 y\ (9) 

Z g = Z\ 9 {Z^I\ (10) 

z 9 = z? 9 {z 3 )-v\zt)-\ (11) 

z g = z^ 9 {z 3 y x i\z 2 )-\ (12) 



2 For simplicity we set the t' Hooft mass fi = 1 in eq. (5) below. 
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Within the MS scheme each RC does not depend on dimensional parameters 
(masses and momenta) and can be represented as follows 

^) = i + £^P, (13) 

n=l C 

where h = g 2 /(167r 2 ) and the parameter e is related to the running space time 
dimension D via D = 4 — 2e. Given a RC Z(h), the corresponding anomalous 
dimension is defined as 

^^a^_ £ !( „ +1 , 

d/i 2 dh ^ 



Z 2 

dimension as the QCD /^-function 



Customarily, one also defines Z h = Z 2 and refers to the corresponding anomalous 



(3(h) = 2 lg (h) = 2h = -E (3 n h^\ (15) 



?1=0 



dz^(h) 

= ^ g \ri) = zn 
Eqs. (9-12) imply that 

(3 = 7i 9 - 2 73 , (16) 
(3 = 2 7l 39 - 3 73 , (17) 
(3 = 2 7 r 9 - 2 73 c - 7 3, (18) 
(3 = 2 7l ^ 9 - 2 72 - 73 . (19) 

As is well-known there is a one-to-one correspondence between an anomalous 
dimension and the corresponding RC. For instance, Zh obeys an equation 

(-e + (3(h)) ^ -(3(h) (20) 



which leads to 

dh (3 

In general case Z depends on \x through both h and £ L and an analog of eq. (20) 
assumes the form: 



hgZh = JT — 



(-e + m)h d ^ + ,s(hn L d -^--m. m 

Eq. (21) can be easily utilized to reconstruct Z from 7 and (3. As anomalous 
dimensions are more compact than corresponding RCs in what follows we will 
write explicitly only the former. 
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3 Calculations and Results 



Relations (9-12) demonstrate that a minimal set of the RCs necessary to recon- 
struct all coefficients of the bare QCD Lagrangian (5) consists of Z m , all three 
wave-function RCs Z3, Z2, Z%, and the coupling constant RC Z g or, instead, at 
least one from the collection ZY, V G {3g, 4g, ccg, ipipg}- Taking into account 
that Z m , Z 2 and Z g are known with four loop accuracy from the works [1-4], one 
is left with just two specific RC to compute, say, Z 3 and Z%. 

At present there are basically two different ways to perform RG calculations 
at the four-loop level. Both approaches make use of the method of Infrared Re- 
arrangement (IRR) Ref. [18] in order to set zero (possibly after a proper Taylor 
expansion) masses and external momenta. Both eventually employ the traditional 
integration by parts method to compute the resulting Feynman integrals 3 . 

The first one, pioneered in the yearly works of Dubna group [24-26], amounts 
to adding an artificial mass or an external momentum to a properly chosen prop- 
agator of a given Feynman diagram before the expansion in masses and true 
external momenta is made. The artificial external momentum has to be intro- 
duced in such a way that all spurious infrared divergences are removed and the 
obtained Feynman integral is calculable. In practical multiloop calculations the 
condition of absence of the infrared divergences leads to unnecessary complica- 
tions and, in some cases, even prevents from reduction to the simplest integrals. 
The problem was solved with elaborating a special technique of subtraction of IR 
divergences — the .Reoperation [27-29]. This technique succeeds in expressing 
the UV counterterm of every (L+l)-loop Feynman integral in terms of divergent 
and finite parts of some L-loop massless propagators. 

In the second approach the infrared rearrangement is performed by introduc- 
ing a single auxiliary mass to all propagators in each Feynman diagram at hand 
[1,30,31]. No IR divergences can ever appear due to absence of any massless 
propagators. Next, after a proper expansion in all the particle masses (except the 
auxiliary one) and external momenta is performed. The resulting integrals are 
completely massive tadpoles, i.e. Feynman integrals without external momenta 
and with only a single mass inserted in all the propagators. 

In our calculation of Z 3 and Z\ we have used the first, "massless" approach. 
It proved also to be more convenient to compute the RC Z 1 ^ 9 instead of Z% and 
then to find Z 3 from eq. (11). 

The four-loop diagrams contributing to the ghost propagator and to the ghost- 
ghost-gluon vertex to order a\ (altogether about 35000) have been generated with 
the program QGRAF [32], then globally rearranged to a product of some three- 
loop p-integrals with a trivial (essentially one-loop) massive Feynman integral 

3 Though the technical implementations could be quite different, cf. works [19-23]. 
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and, finally, computed with the program MINCER [33,34]. The total amount 
of CPU time needed to compute RC and Z c ^ g was about a month of work 
of a standard PC with an Athlon XP 2000+ processor. For testing purposes we 
have also computed the RC Zf q9 , which has required an almost double amount 
of calculational time 4 . 

Our results for the anomalous dimensions jl° 9 and 73 read 



(22) 



(23) 



(7f*) S 



5427 _ 1053 _ 2 135 _o 135 



64 



64 



32 



16 



n f , 



(24) 



(7? 



cc 9\ 



+ 
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384 ' 32 
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64 
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128 



, 243 3 



+ n / 



128 
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945 



64 ' ^ 128 " 128 
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~2gg-eL--g-Cs&- T C4& 



64 Cs£l + 12 g Cs^Z 



453 
132" 
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54 



(25) 



(73)o = -^ + ^L, 



(26) 



285 
76" 



9 , 5 
16^+4^' 



(27) 



(7 3 c ) 2 = 



15817 243 



32 



459 
~32 



'A 



+ n f 



64 

32 S3 ^ L 64 
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32 



CI 



24 



16 



35 2 
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36 p 



(28) 



4 These figures should be understood as effective ones; that is in reality we have used a sort 
of trivial parallclization by distributing diagrams between a few PC's. 
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(73 C )s = 
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+ 
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Finally, a use of eq. (18) immediately leads us to 



(29) 



13 3 2 
(73)0 = -y + 2^+^/3, 



531 99 , 9 , 61 
(73)i = -— + y^ + i £ L +yn /) 



(30) 
(31) 
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For an important particular case of the Landau gauge £ L = we get: 



13 2 
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(73)2 
(73)3 = 
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or, numerically, 

7s = 



= h(Q.5 - 0.666667/1/) + /i 2 (66.375 - 10.166771/) 
+ /i 3 (915.963 - 186.86^ + 7.96296^) 
+ /i 4 (19920.2 - 4692.27n/ + 331.285nJ - 0.169134nJ). 



(33) 

(34) 
(35) 



(36) 



(37) 



4 Applications 

In this section we consider some applications of our results. The case of the 
massless QCD with the Landau gauge fixing is understood in both subsections. 



4.1 Scheme- invariant Gluon and Ghost Propagators in 
NNNLO 

In general case a (multiplicatively renormalizable) Green function G depends on 
both a renormalization prescription (scheme) and the choice of the normalization 
scale fi. In many cases it is more convenient to deal with the scheme and scale 
invariant version of G which we will denote as G. Given the RG equation for G 

fi 2 -^G(h, n) = (V A + (3(h)h^r) G = 7 (/i) G(h, p), (38) 
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then a formal solution for G reads 



G = G(h,Mf(h), /(^expi/^lj}, 



(39) 



f(h) = (/i)*{l + (7i-A7b)fc 

+ ^ [(71 - Ato) 2 + 72 + /9i 2 7o - A71 - 



h 2 



+ 



11 2 

-(fi - Pilof + ^(71 - /5i7o)(72 + Pi 7o - A71 - /3 2 7o) (40) 
+ ^ ( 73 - A 3 7o + 2/3i/3 2 7o - &7o + /3?7i - &7i - A72)] ^ + ^(^)} • 



Here 7$ = 7i//3o, A — A/Ah (i=l,2,3) , the coefficients 7,, A are defined in 
eqs. (14,15). 

Now, after writing the gluon and ghost propagators in the form 



sab 



D ab (q) = 



Sab 

D(-q 2 ), A ab (q) = ^ A(-q 2 ) (41) 
-q 2 



and a use of explicit expressions for the very propagators from [35] we arrive at 
the following NNNLO predictions for the asymptotic behavior of the scheme and 
scale invariant functions D and A at large Euclidean Q 2 = —q 2 in the MS scheme. 
First, for the case of pure gluodynamics (rif = 0) 



h% D~\-q 2 



|n/=0 



n 25085 , , 2 
1 — ^r^h+ h 2 



2904 

, , r 141629801206331 80968605 . 

I h + Cs 

L 16326706176 42592 s 

= 1- 8.63809/1-186.819 /i 2 -5686.55 h 3 , 



412485993 9747 _ ' 
+ —(3 



1874048 
477315 ' 

704 ^. 



352 



(42) 



h& A~\-q 2 ) 



|n/=0 



n 5271 , , 2 

1 h+ h 2 

1936 



615512003 5697 
7496192 + 704 ^ 3 



+ h 



430343889400537 674654895 , 73845 , 



130613649408 1362944 
1-2.72262 /i-72.3825 /i 2 -2482.24 h 3 . 



352 



(43) 



In order to illustrate the rif dependence we give below the results for rif = 3 
and rif = 6 (to save space only in the numerical form) 

hh D'\-q 2 )\ nf=3 = 1 - 5. 18056/i - 85.0853 h 2 - 2178.1 h 3 , (44) 



hA A (-q )\ n =3 = 1 - 2.78472/i - 52.591 /i 2 - 1359.11/^, 



(45) 
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h& b~\-<?)\ nt =fi = 1 - 0.857993 h + 16.6153 h 2 + 220.455 h 3 , (46) 

}Jg K-\-q 2 ) | n/=6 = 1 - 3.27934 h - 31.5908 h 2 - 372.071 h 3 g. (47) 

Note that in eqs. (42-47) the coupling constant h should be understood as 
h^ 2 = -q 2 ). 

4.2 Four-Loop Anomalous Dimension of the Composite 
Operator A 2 

Recently there has been a lot of activity in studying the possibility of a condensate 
in Yang-Mills-theory of mass dimension two (see, e.g. recent works [36,37] and 
references therein). The relevance of the operator A 2 = A® A® in the Landau gauge 
in that context has been widely discussed. In this connection a thorough inves- 
tigation of the renormalization properties of the composite operator A 2 = A^A^ 
has been carried out in [38,39]. In particular, the author of [38] has discovered 
by explicit three loop computation a remarkable relation 5 

-2 7A2 |^ =0 = /3-7 3 (48) 

expressing the anomalous dimension 7^42 of the operator A 2 in terms of the f3- 
function and the gluon field anomalous dimension 73. An all orders proof of (48) 
have been later constructed in [40] using algebraic renormalization methods. As 
both ingredients of (48) are known now at four loops one arrives at 
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+ 



+ 



h 3 
h A 



h [35 

■75607 
.~64 

■29764511 
. 1536 

+ n f (- 



2 

243 
l32~ 

99639 



+ h 2 
Cs + n f ( 
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162 
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(49) 



5 Discussion and Conclusion 

Calculation of the two missing RCs Z 3 and completes renormalization of the 
QCD Lagrangian at four loops. An important issue relevant for any calculation 

5 We have adjusted the coefficients in (48) to our notations. 
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of such complexity is the correctness of the obtained anomalous dimensions. The 
following comments are in order. 

• The FORM program MINCER used by us to compute three loop massless 
propagators was developed more than a decade ago and has been since 
heavily cross-checked in a number of various multiloop calculations. 

• At three loop level we have full agreement with the results of [4,5]. This 
checks our way to use -Reoperation because the three loop results of [4, 5] 
have been obtained with direct application of MINCER to three-loop prop- 
agators. In the present work the three loop contributions come exclusively 
from two loop propagators. 

• The leading rif behaviour of the quark, gluon and ghost anomalous dimen- 
sions as well as the anomalous dimensions of the quark-quark-gluon and 
ghost-ghost-gluon vertices was investigated in all orders of perturbation 
theory in the work [41]. At the four loop level the predictions of [41] are in 
full agreement to the corresponding leading rif pieces of our results 6 . 

• We have also computed the 0(aj) anomalous dimension of the quark-quark- 
gluon vertex and found that it satisfies eq. (19) as it should. 

• As is known from [42, 43] the ghost-ghost-gluon vertex is unrenormalized in 
the Landau gauge, that is 

7rk=o = 0. (50) 
Eq. (50) is in obvious agreement to eqs. (22 - 25). 

• We have not performed a direct calculation of the gluon wave function RC Z 3 
but rather extract the result from the /9-function of [1] and the anomalous 
dimensions 73 and j^° 9 . Thus, an independent reevaluation of either 73 
or/and the /3-function at four loops is highly desirable. 

An interesting feature of the scheme and scale invariant functions D and A 
is the full absence of the irrational constant (4 as illustrated in eqs. (42,43). It 
stems from a noni-trivial mutual cancellation of the terms proportional to £4 which 
enter into the ingredients — the Green function and its anomalous dimension — of 
the definition of G (see eq. (39)). We have checked that this is true also for the 
SU(N) color group and a generic value of rif. 

Finally, a comment about the gauge group dependence of our results. Pre- 
ferring compact and readable formulas to huge and general ones, the author has 
deliberately formatted all results in this paper for the most practically important 

6 We thank John Gracey for this comment. 
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case of the SU(3) color group. In reality all calculations have been carried out 
for a little bit more general case of the SU(N) color group. Full expressions of 
RCs (and the corresponding anomalous dimensions) describing the renormaliza- 
tion of the Lagrangian (5) in the general covariant gauge and for the SU(N) color 
group are available (in a computer-readable form) in http:/ /www-ttp.physik.uni- 
karlsruhe.de/Progdata/ttp04 /ttp04-08 / . 

Acknowledgments: I would like to thank Michael Czakon and York Schroder 
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want to thank Mboyo Esole for information about the work [37]. 
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Note 1 

Just before the submission of the manuscript for publication we have been 
informed that our main assumption — the validity of the result for four-loop QCD 
beta-function first obtained in [1] - is confirmed by a completely independent 
calculation [44]. 
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